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Abstract 

We investigate J /if) production at a Photon Linear Collider (PLC) in the framework of 
NRQCD factorization approach. It was shown that color octet contribution is dominating 
in charmonium production at PLC. The main fraction of J ftp comes from the single and 
double resolved photon processes. Expected cross sections of charmonium states production 
are large enough even at px > 5 GeV to expect several hundreds of J/ip's through the 
leptonic decays channel for a projected integral PLC luminosity of 50/6 -1 
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1 Introduction 

The nonrelativistic QCD (NRQCD) factorization approach, developed in the last few years, is the 
most successful method for the analysis of heavy quarkonium production and decay processes [jlj . 
In the framework of this formalism cross section and decay rate of the process involving heavy 
quarkonium is factorized into the short distance part (coefficient functions) and the long distance 
one. For the cross section a(A + B — > H), where H stands for quarkonium, the formalism states: 

a ( A + B H) = 5>(4 + B -> QQ[n})(0\O H [n}\0). (1) 

n 

In (1) a(A + B — > QQ[n}) is the cross section of the production of the heavy quark-antiquark 
pair state [n]. This part of cross section can be calculated perturbatively. The matrix element 
(0\O [n]|0) stands for the evolution of the state [n] into the particular hadronic state H. It is 
defined by long distance dynamics and cannot be computed perturbatively. The long distance 
matrix elements scale as a power of v - the relative velocity of the heavy constituents in the 
quarkonium. So, in the NRQCD factorization approach (FA) the cross section of heavy quarko- 
nium production is considered as a double series in the QCD coupling constant and the relative 
velocity v. The relative importance of long distance matrix elements in powers of v can be 
estimated using the NRQCD velocity scaling rules |fj]. 

In the NRQCD factorization approach the quark-antiquark color octet intermediate states 
are allowed to contribute in heavy quarkonium production and decay processes in higher order. 
By such a way the Color Octet Mechanism (COM) takes into account the complete structure 
of the quarkonium Fock space while in the color singlet model (CSM) only the dominant Fock 
state is considered @. The shape of the px distribution of the 3 S[ 8 ^ octet state production cross 
section indicates that J / ip and ip' production at large px observed at the Tevatron (FNAL) can 
be explained in the framework of FA ||. |5| . 

Despite the obvious successes of the COM some problems remain unsolved. In particular, the 
theoretical predictions disagree with the J/ip and ip' polarization data at fixed target energies 
H 0] and the COM estimation for the yield ratio of x& an d Xc2 states remains too low ||. These 
discrepances indicate that higher twist corrections might contribute significantly in production 
of charmonium states at low pt and should be considered equally with the COM contributions 
||. The latter underestimate the J /if) photoproduction cross section at large values of z in HI 
and ZEUS experiments at HERA (z = Ej/^/E^ in the proton rest frame)- the most problematic 
issue of the COM @. 

The NRQCD factorization approach has predictive power when the relative velocity v is small 
enough, allowing to describe with reasonable accuracy several experimental data by taking into 
account only a few lower order terms in powers of v. Hence it becomes necessary to check the 
universality of the FA, stating that the values of long distance matrix elements, extracted from 
the experimental data for different processes, should be the same. Unfortunately, due to the 
rather large theoretical uncertainties at the present time, the existing experimental data does 
not allow to extract the values of the long distance matrix elements with enough accuracy. 

This has motivated us to look for processes with less theoretical uncertainties to test the color 
octet mechanism. The measurements of the J / ip polarization in unpolarized hadron-hadron col- 
lisions and electroproduction [§, [U], [11], [12| can be used for these purposes as well as the asym- 
metries in J ftp production JT2J. Unlike hadron-hadron collisions, existing J/ip photoproduction 
data can be explained within the CSM taking into account higher order QCD corrections fll3| . 
Moreover, color octet contributions overestimate the experimental data drastically |J. In this 
sense, the study of J ftp production processes in non-hadron collisions can clarify the existing 
contradictory picture. 
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The collision of high-energy, high-intensity photon beams at a photon linear collider, obtained 
via Compton backscattering of laser beams off linac electron beams, should provide another 
possibility to investigate the heavy quark-antiquark bound states production processes. Based 
on a e + e~ linear collider, the PLC will have almost the same energy and a luminosity, i.e. a cm. 
energy of 500-1000 GeV and a luminosity of the order of 10 33 cm" 2 s _1 ||T4| . 



In the present paper we study J/ip production at a photon linear collider. The color octet 
contribution in J/ip production at PLC was calculated recently fT5| , where only the direct 
(unresolved) photon-photon collisions were considered. We will show that the resolved photon 
processes are more significant because there exist diagrams in lowest order of perturbative QCD 
leading to the subprocess cross section, non decreasing when the invariant mass of photon- 
photon system increases. In the next section we present all possible subprocesses in lowest 
order in perturbative QCD and fix the color octet long distance matrix elements needed. The 
numerical results are given in Sec. III. It is shown the that color octet contribution in resolved 
photon subprocesses is essential and the expected cross section of charmonium production is 
large enough even at high (> 5 GeV) allowing to extract valuable information about color 
octet long distance matrix elements. 



2 J/ip Production Subprocesses and COM Parameters 

For the calculation of the J/ip production cross section at PLC we consider the following sub- 
processes: direct 77 collision, single and double resolved photon processes, when the heavy 
quark-antiquark pair is produced via hadronic components of the photon, which are described 
in terms of the photon structure functions fl6| . 

To lowest order in a em and a$ the following subprocesses contribute: 

7 + 7 -> (cc) + -rig) (2) 
7 + 7(0,9) (cc) +9,Q (3) 

7(9, q) +i(g,q) -> (cc) +g,q (4) 

where j(g, q) denotes a gluon or quark component of the photon, (cc) stands for any possible 
quark-antiquark state. 

The cross section of J/ip production in the NRQCD factorization scheme can be presented 

as 



<j m = <r( J/ip) dir + E Br(xcj->J/i>X)a XcJ + Br(il/->J/il>X)a i , (5) 

J=0,l,2 

Here Br((cc) — ► J/ipX) denotes the branching ratio of the corresponding (cc) state into J/ip. To 
the production of each state of quarkonium both color octet and color singlet states contribute, 
as in the case of direct J/ip production 

a(J/^) dir = a-f + a* m = a(J/^ 1 + J2^(QQ[ 2s+1 L S I })(0\O J ^C s+1 L J )\0) (6) 

where the sum runs over the states 3 Pq±2^ an d 3 Sf. In the framework of the NRQCD FA 
velocity expansion we consider only the contribution of the dominant sets of color octet states 
in the direct production of S and P states of charmonium. 
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In direct 77 interaction in lowest order only the two states of heavy quark-antiquark pairs 
can be produced : 3 Si color singlet (J/ip) and 3 Si color octet For the differential cross 

section of the singlet state production we have: 

at 6 6 

s 2 (s - M 2 ,,,) 2 + i\i - M 2 ^) 2 + u 2 (u- M 2 ^) 2 

s\s-M^f{i-M^f{t-M^Y [ } 

where s, t and u are the standard Mandelstam variables. 

The differential cross section for the J /if) production through color octet 3 Si state can be 
obtained from (0) by the replacement Q 

(O j ^( 3 Si)) - ±£L(0>'+(> Sl )) 

We calculated all possible subprocesses for single resolved photon reaction. Our results are 
in agreement with ||. The subprocess for (cc) states production in double resolved photon case 
are same as for J/ip production in hadron-hadron collision and were calculated by Cho and 
Leibovich ||. 

All color singlet long distance matrix elements are related to the radial quarkonium wave 
functions at the origin and their derivatives. We have used the Buchmiiller-Tye wave functions 
at the origin, tabulated in |1. The number of color octet long distance matrix elements can be 
reduced by using the NRQCD spin symmetry relations: 

(0|Of( 8 Pj)|0> = (2J+l)(0|Of( 3 P )|0), (8) 
(0|O^( 3 ^)|0> = (2J+l)(0|O 8 ^( 3 5 1 )|0). (9) 

These relations are valid up to order v 2 . 

After using these relations only three independent matrix elements remain - (Og ( 3 jSi)), 
(Cg^( 3 Po)) and (Og^ ( x S )} , which give the main contribution to the direct J/ip production 
cross section. The values of these parameters are extracted from the data of direct J/i[) produc- 
tion at CDF [10[| . We used the results of paper HlOj , where a more detailed numerical analysis 



was carried out: 

) J /^3(v\|n\ _ -I nfi 4- n 1 /I +1-05 



( Si)|0) = 1.06 ± 0.14™ • 10" 2 GW 3 , (10) 



(O|0 8 J/ YSo)|O> + ^(0|0 8 J/ ^( 3 Po)|0) = 3.9 ± 1.151J-J? ■ !Q- 2 GeV 3 (If 

mi 



Further we will assume that (0|C 8 7/,/ '( 1 5o)|0) = 3.5/m 2 c (0\O^( 3 P )\0) . For the indirect J /if) 
production via the Xcj states decays we use the following value fitted from the CDF data H: 



(0$ cl ( 3 Si)} = 9.8 • 10- 3 GeV 3 . (12) 

We would like to to mention that all these values contain large uncertainties (for more detailed 
analyses see || [TO], [L7|]) and we use them only to estimate the expected cross sections of the 



production of different charmonium states at PLC. 

2 The additional factor 4/3 coming from different color factors of considered processes is omited in the paper 
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3 Numerical Results 



At PLC the energy distributions of the colliding photon beams are not monochromatic. One of 
the main characteristics of the Compton scattering process is the dimensionless variable x, 



x 



4E b u 



mi 



cos 0/2. 



Here £7& is the linac beam energy, uq is the energy of the laser photons and 9 is the angle between 
the directions of the electron and laser beams. The maximum energy of the obtained photon 
beam depends on x, 



E h 



x 



1 + x' 

and increases with increasing x. The value x = 4.83 is optimum for PLC to prevent e + e~ pair 
creation in a collision of laser and backscattered photon beams f|XQ|l . In our calculations we use 
x = 4.8. This means that the maximum energy fractio carried by backscattered photon, ut/Eb, 
is about 0.83. 

The luminosity distribution of the photon beam depends sensitively on the conversion dis- 
tance due to non-zero angles of photon-electron beam scattering. The conversion distance is 
the distance from the laser and electron beam intersection point to the final interaction point. 
The luminosity distribution depends also on the electron beam size and shape. Assuming round 
Gaussian linac beams, the luminosity spectrum is characterized by a geometrical factor p, the 
ratio of the intrinsic transverse spread of the photon beam to that of the original one : 



3 . 61v ^T(±)(Ji_rv 

V cm / \1 ev J \ nm 



-i 



(13) 
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Figure 1: Spectral luminosity of 77 collisions; the 
dashed line corresponds to the value p = 0, the solid 
line to p = 0.6 and the dashed-dotted line p = 2.4. 



here / is the conversion distance and #0 is 
defined by the electron beam energy, 9 = 
m e \/l + x/Eij. In fig.l the spectral lumi- 
nosities for unpolarized photon and electron 
beams are presented for different values of 
p are presented. The parameter k is the ra- 
tio of the number of electrons scattered on 
the laser beam to the number of electrons in 
each electron bunch . With p increasing the 
spectral density becomes narrower and more 
energetic photons dominate, meanwhile the 
total photon photon luminosity is decreas- 
ing. To achieve a maximum J / ip production 
cross section the geometrical factor p should 
be as small as possible. However, in this case 
the transverse separation of the electron and 
the secondary photon beams becomes im- 
possible and large background is expected 
from the process —>■ eZ —>■ ebb with the 
subsequent decays of b(b) quarks into J/ip. 
Hence we use p = 0.6 which is small enough 
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to achieve an observable J/ip production cross section (see below) as well as large enough to 
eliminate electron-photon collision background [ED . 



Table 1. J/ip production cross sections in direct 77 interaction at PLC for y/s e + e - = 500 
GeV and 800 GeV. The cross sections are calculated for two different values of the geometrical 
factor p =0 and 0.6. 



p 


^=500 GeV 


v^=800 GeV 


singlet 


octet 


singlet 


octet 


0. 


2.79 ■ 10- 5 n6 


2.71 ■ 10- 5 nb 


1.21 • 10- 5 nb 


1.18 • 10~ 5 nb 


0.6 


1.74 • 10~ 7 nb 


1.69 • 10~ 7 n& 


3.10 • 10- 8 nb 


3.00 • 10~ 8 nb 



In Table I J/ip production cross sections for two different values of the geometric factor, 
p = and 0.6, are presented. We have used the color octet long distance parameters presented 
in the previous section. As one can see from table I, the total cross section of J/ip production 
is decreasing with increasing e + e~ collision energy. The reason for this is that the production of 
the 3 Si state (singlet or octet ) in direct 77 interaction falls off as 1/s 2 with increasing photon- 
photon invariant mass. For p = the luminosity of low energetic photons is higher than for large 
values of p and correspondingly the J / ip production cross section is higher. The J / ip production 
cross section in direct photon-photon collision for p = 0.6 is practically unobservable even at a 
high projected e + e~ collider luminosity of 100/6 -1 . 

As for the resolved photon subprocess, unlike direct 77 collisions, there exist diagrams with 
only vector particles (gluon) in the t-channel in the lowest order of perturbative QCD. These 
diagrams result in a cross section of (cc) pair production which is not decreasing with increasing 




1.5 2 2.5 3 3.5 4 4.5 5 1.5 2 2.5 3 3.5 4 4.5 5 

Pt ( GeV ) p T (GeV) 

Figure 2: The J ftp production cross sections at ^Js e + e - = 500 GeV, a) in single resolved photon 
processes, b) in double resolved photon processes. 

In fig. 2 we present the direct J/ip production cross sections through different (cc)-pair states 
for the single resolved (fig. 2a) and double resolved (fig. 2b) photon subprocesses. We have used 
m c — 1.5 GeV for the charm quark mass and GRV pi] photon distribution functions evaluated 
at the factorization scale Q 2 = \Jpx + 4m^. The QCD coupling constant as is normalized at 
the same scale. We present the numerical results for the cross sections of ISq and 3 Pj states 
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production for px > 1.5 GeV in order to avoid the consideration of mass singularities and to 
deal with experimentally observable events. The J/tp production cross sections in both single 
and double resolved photon processes are higher than in direct 77 collisions as can be seen when 
comparing fig. 2 and table I. The color octet contribution is dominating in J/ip production in the 
resolved photon case. At small values of px < 5 the color octet states 1 5'o and 3 Pj give main 
the contribution in direct J/ip production. As in hadroproduction, the 3 S[ 8 ^ octet state becomes 
essential only at large values of px (> 5 GeV) in the double resolved photon processes, as can 
be seen from fig. 2b. 

Note that a nondecreasing cross section of J / ip production, with increasing 77 invariant mass, 
can be obtained also in direct (unresolved) photon-photon reactions when considering higher 



order diagrams in pQCD with two-gluon exchange in the t-channel |[22|| . But the contribution 
to the total cross section of J/ip production from such processes, cr to t{ll J/^X) ~ 2 • 10~ 2 nb 
is smaller than color octet contribution in single and double resolved photon subprocesses, 



and decreases rapidly at large transverse momenta. 

In fig. 3 we present the direct J /if) production total cross sections in dependence on the cut 
off value in px, for different angular cuts in laboratory frame (e + e~ c.m.s. system). Fig. 3a 
corresponds to ^Js e + e - = 500 GeV and fig. 3b to 800 GeV. The leading contribution in the cross 
sections with cuts cos# < cos 30° or cos 45° comes from the double resolved photon processes, 
the contribution of single resolved photon subprocesses is negligible. In single resolved photon 
subprocesses (cc) octet states are produced with large longitudinal momenta. As in the case 
of photon-proton collisions, color octet states are mainly produced at the kinematical endpoint 
z — 1, as in the case of J/ip photoproduction [[J; here z = (p l2 ■ Pj/^)/{p l2 ' Pn)-> where p 7l is 
unresolved photon momentum and p 12 is the momentum of the photon which splits into a flux of 
quarks and gluons. This means that the J/tp carries almost the whole energy of the unresolved 
energetic photon and is produced at small angles relative to the beam direction. 




1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

p T min (GeV) P T min (GeV) 



Figure 3: The J/tp production cross section versus minimal px cut for different angular cuts. The solid 
line corresponds to the total cross section, the dashed line to cos 9 < cos 30°, and the dashed-dotted 
line to cos6» < cos 45°; a) yjs e+e - = 500 GeV and b) yjs e+e - = 800 GeV. 

The double resolved photon process is analogous to J/ip hadroproduction. At large px 
the contribution from the 3 S[ 8 ^ color octet state dominates. Table II presents the integrated 
total cross sections and 3 S' 1 8 ^ channel contribution for various p™ m cuts. All cross sections are 
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calculated at an angular cut off cos 8 < cos 45°. The cross sections are large enough to observe 
several hundred J/ip events produced at px > 7 GeV through leptonic decay channels for the 
projected integral PLC luminosity of £ = 50 fb^ 1 . This would allow to extract the long distance 
matrix element (0$^( s Si)) with sufficient accuracy []. 

Table II. The J/ifj production cross sections (total and through 3 S[ 8 ^ state), cos 6 < cos 45°. 



Pt 


5 GeV 


6 GeV 


7 GeV 


total 


6.2 • 10~ 4 n6 


3.2 • 10~ 4 n6 


1.8 • 10- A nb 


3g(8) 


3.8 ■ 10- 4 nb 


2.2 ■ 10~ 4 nb 


1.4- lfr 4 n6 



Fig.4 presents the J/ip production cross section resulting from the radiative decays of the Xcj 
state. The cross sections are calculated with an angular cut cosO < cos30°. At subleading order 
in the relative velocity expansion, 0(v 5 ), only one color octet state, 3 S[ 8 \ contributes to the pro- 
duction of Xcj (the corresponding color octet long distance parameter was extracted only from 
the Tevatron data The next corrections are already of order 0(v 9 ) H). In fig.4 the color octet 
(dashed-dotted line) and color singlet (dotted line) contributions to Xcj meson production are 
presented. 
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1.5 



2.5 



3.5 



4.5 
p T min 



5 5.5 

(GeV) 



°(p T >p T min , ^S=500 GeV) (nb) 

4 Conclusions 

In the present paper we consider J / ip meson pro- 
duction at PLC. We calculated cross sections for 
the production of different color octet and color 
singlet heavy quark-antiquark states in direct and 
resolved photon processes. In the resolved pho- 
ton subprocesses already at the lowest order of 
perturbative QCD the diagrams arise which lead 
to the non decreasing production cross sections 
for some quark-antiquark states with increasing 
photon-photon c.m.s. energy. 

The contribution of color octet states to J/ip 
production is dominating. It was shown that at 
large pt (pt > 5 GeV) the main contribution to 
J/ip production comes from the 3 S± state. The 
cross section is even at such values of transverse 
momentum large enough to observe several hun- 
dred J/ip's through the leptonic decay channel at 
a projected PLC luminocity 50/6" 1 . This should 
allow to extract the value of the corresponding 
matrix element. The cross section at smaller val- 
ues of transverse momenta can also give valuable 
information about the color octet long distance matrix elements. The production cross section 
of J/ip mesons through Xcj states was also calculated. It is worth mentioning that background 



Figure 4: Cross section for J/ip production re- 
sulting from radiative Xcj decay versus mini- 
mal pt cut. The dashed curve corresponds to 
the color singlet cross section, the dashed-dotted 
line shows the 3 iSj 8 ^ contribution, the solid curve 
represents their sum. 



5 We mean only statistical significance and do not take into account theoretical uncertainties 
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from fe-quark decays are smaller or at same order as direct J / ip production p3| but the analysis 



of possible backgrounds and ways of their rejection are subject of an additional study. 
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